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Single crystals of the fiUed-skutterudite compound NdOs4Sbi2 have been investigated by means 
of electrical resistivity, magnetization, and specific heat measurements. The NdOs4Sbi2 crystals 
have the LaFe4Pi2-type cubic structure with a lattice parameter of 9.3 A. Possible heavy-fermion 
behavior is inferred from specific heat measurements, which reveal a large electronic specific heat 
coefficient 7 ~ 520 mJ/mol-K^, corresponding to an effective mass m* « 98 rUe. Features related 
to a ferromagnetic transition at ~ 0.9 K can be observed in electrical resistivity, magnetization and 
specific heat. Conventional Arrott-plot analysis indicates that NdOs4Sbi2 conforms to mean-field 
ferromagnet ism . 



PACS numbers: 75.40.Cx, 71.27.+a, 71.70.Jp, 71.70.Ch 



I. INTRODUCTION 

The filled skutterudite compounds have the chemical 
formula MT4X12, where M is an alkali (Na or K), al- 
kaline earth (Ca, Sr, Ba), lanthanide or actinide atom; 
T is a transition- metal atom (Fe, Ru, Os); and X is a 
pnictogen atom (P, As, Sb). The compounds crystal- 
lize in a LaFe4Pi2-type structure with the space group 
Im3.^ Due to the strong hybridization between f- and con- 
duction electrons and their unique crystal structure, the 
lanthanide- and actinide-based filled skutterudite mate- 
rials display a wide range of strongly-correlated-electron 
phenomena, such as BCS-like superconductivity (e.g., 
PrRu4Sbi2)ri^ heavy fermion behavior (e.g., PrFe4Pi2)(^ 
heavy fermion superconductivity (e.g., PrOs4Sbi2)r*^ 
ferromagnetism (e.g., PrFe4Sbi2)f2, metal-insulator tran- 
sitions (e.g., PrRu4Pi2))^Kondo-insulator behavior (e.g., 
UFe4Pi2 and CeFe4Pi2)f^ valence fluctuation behavior 
(e.g., YbFe4Sbi2)r2*ii and non-Fermi-liquid behavior 
(e.g., CeRu4Sbi2)i2^ 

Previous studies have shown ferromagnetic order- 
ing in Nd-based filled skutterudites. Measurements of 
NdFe4Sbi2 revealed ferromagnetic order below 16.5 K 
with a Nd ordered moment of 2.04 /Lte and a Fe 
collinear moment of 0.27 /xb"^ Lower ferromagnetic 
transition temperatures were found in NdFe4Pi2 at 
2K^ NdRu4Pi2 at 1.6 K,!^ NdRu4Sbi2 at 1.3K,2-3 and 
NdOs4Sbi2 at 0.8 K.^^ Since NdOs4Sbi2 displays ferro- 
magnetism with a low Curie temperature and its neigh- 
boring compound PrOs4Sbi2 shows heavy fermion be- 
havior and unconventional superconductivity)SiiLiS4i2i2£ 
possibly involving triplet spin pairing of electrons, there 



is a strong likelihood that PrOs4Sbi2 is near a ferromag- 
netic quantum critical point. Thus, by thoroughly char- 
acterizing the physical properties of NdOs4Sbi2, a deeper 
insight into the unusual behavior of PrOs4Sbi2 may be 
attained. Earlier studies of the compound NdOs4Sbi2 
only reported the results of structural refinements^ and 
the value of the Curie temperature^ In this report, we 
present a new and detailed investigation of NdOs4Sbi2 
single crystals, including measurements of X-ray diffrac- 
tion, electrical resistivity, magnetization, and specific 
heat. We also discuss possible heavy fermion behavior 
in this Nd-based filled skutterudite compound. 



II. EXPERIMENTAL DETAILS 

NdOs4Sbi2 single crystals were grown in a molten Sb 
fiux as described previously,'^'^^-^ using high purity Nd, Os 
(3.5N), and Sb (6N). X-ray powder diffraction measure- 
ments were performed with a Rigaku D/MAX B x-ray 
machine on a powder prepared by grinding several single 
crystals, which indicated single phase NdOs4Sbi2 with a 
minor impurity peak of Sb (< 10%). The crystals had 
a LaFe4Pi2-type BCC structurefi with a lattice param- 
eter a — 9.30 A. Two single crystals of similar dimension 
were selected for single crystal X-ray diffraction measure- 
ments. The data were collected on a four-circle Nonius 
Kappa diffractometer at 296 K using Mo radiation 
(A = 0.071073 nm). No absorption corrections were nec- 
essary because of the rather regular crystal shape and 
small dimensions of the investigated crystals. The struc- 
ture was refined with the SHELXS-97 programi2^ 
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Electrical resistivity p(T, H) was measured using the 
standard 4-wire technique in a Quantum Design PPMS 
system and in a ^He-^He dilution refrigerator in fields 
up to 8T. The low temperature (0.02 K - 2.6 K) and 
high-field (8 T - 18 T) p{T, H) measurements were per- 
formed in the National High Magnetic Field Laboratory 
at Los Alamos National Laboratory. The electrical cur- 
rent applied to the sample was perpendicular to the ap- 
plied magnetic field, which was along the [001] direction, 
in all p(T, H) measurements. Measurements of p(T, P) 
were made under hydrostatic pressures up to 28kbar in 
a beryllium-copper piston-cylinder clamp^^ and a ^He 
cryostat. The pressure was determined inductively from 
the pressure-dependent superconducting transition of a 
Pb manometer. 

DC magnetic susceptibility from 2 K to 300 K was mea- 
sured in a Quantum Design MPMS SQUID magnetome- 
ter. The magnetization M{H, T) measurements were car- 
ried out in a '^He Faraday magnetometer with a gradient 
field of ^ 0.05 - 0.1 T/cm in external fields up to 5.5 T and 
at temperatures between 0.4 K and 2K. For the Faraday 
magnetometer measurements, several single crystals (to- 
tal mass of 21.3 mg) were combined in a mosaic fashion 
and measurements were performed with magnetic field 
applied along the [001] axis. Specific heat C[T) of multi- 
ple single crystals (total mass of 42.15 mg) was measured 
between 0.5 K and 70 K in a He calorimeter using a semi- 
adiabatic heat-pulse technique. 
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FIG. 1: (a) DC magnetic susceptibility Xdc vs temperature 
T from OK to 10 K measured at 500 Oe. Fits of Xdc(r) to 
a CEF model in which the ground state is either Fg (dashed 
line) or Fg'^' (solid line). Below 10 K, these two lines overlap, 
(b) Xdc^ vs T from K to 300 K. Fits of Xdc(r) and Xdc(T) to 
a CEF model in which the ground state is either Fe (dashed 

-,(2) 



line) or Fg (solid line). 



III. RESULTS AND DISCUSSION 
A. Single crystal structural refinement 

Structural refinement was performed on X-ray diffrac- 
tion data collected from single crystals of NdOs4Sbi2; 
the results are listed in Table ^ The thermal displace- 
ment parameters Un of the Nd atoms are isotropic and 
have large values compared with the Ua for the Os and 
Sb atoms, a common feature in the filled skutterudites. 
The Nd sites are fully occupied, which is not always the 
case for filled skutterudites, such as Pro.73Fe4Sbi2 and 
Euo.95Fe4Sbi2.'''^'' If the NdOs4Sbi2 crystal is considered 
to be a simple Debye solid with the Nd atoms behaving 
like Einstein oscillators, the thermal displacement and 
the Einstein temperature 0e are related by 

U = coth ( ^\ , (1) 

where m^^ is the atomic mass of Nd. For NdOs4Sbi2, 
0E is estimated as ^ 45 K, which is close to the values 
found for thallium-filled antimony skutterudites such as 
Tlo.22Co4Sbi2, Tlo.5Co3.5Feo.5Sbi2, Tlo.8Co3FeSbi2, and 
Tlo.8Co4SbiiSni^ 



B. Magnetic Properties 

The dc magnetic susceptibility Xdc{T) of NdOs4Sbi2 
was measured at 500 Oe and is displayed in Fig. ^a). 
The XdcC^) data (Fig. Ub)) exhibit different slopes at 
high and low temperatures. The linear slope of XdcC^) 
between 65 K and 300 K yields a Curie constant (Ccw 
= (NAAicf[)/(3fcB)) ~ 1.85cm^K/mol, a negative Curie- 
Wiess temperature (Ocw) ^ -43 K, and a effective mo- 
ment HeS ~ 3.84 ^B, which is close to the Nd^+ free- 
ion value of 3.62 /iB. The Curie- Weiss fit to the low- 
temperature range of XdcC^) (Fig- 121^)) gives Ccw 
0.69cm^K/mol, a positive 8cw IK, and a value of 
/ioff ~ 2.35 ^B- The curvature in XdcC^) ^'^ 
infiuence of the crystalline electric field (CEF), and the 
positive Ocw from the low-temperature fit indicates fer- 
romagnetic order developing below IK. 

Although the crystal structure of NdOs4Sbi2 has tetra- 
hedral symmetry (Th)r^ this is only a slight deviation 
from cubic symmetry (0;i). Thus, to simplify the CEF 
analysis, only O^ symmetry was considered. In an ionic 
(localized) model with cubic symmetry, the Nd'^+ ten- 
fold degenerate 3 — 9/2 Hund's rule ground state multi- 
plet splits into a Fg doublet and two Fg (Fg^^ , Fg^') quar- 
tet states. In the treatment of Lea, Leask, and Wolf^^ 
these energy levels and their corresponding wave func- 
tions in cubic Oh symmetry can be parameterized by the 
variables xllw and W, where xllw is the ratio of the 
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TABLE I: Single crystal structural data measured at T = 296 K for NdOs4Sbi2. The crystal structure is LaFe4Pi2-type with 
the space group Im3 (No. 204). The range of X-ray scattering angle is 2° <29 < 80°. 



NdOs4Sbi2 


Crystal size 


64 X 78 X 84 /im'' 


Lattice parameter a [A] 9.3075(2) 


Density p [g/cm'*] 


9.745 


Reflection in refinements 


473 < 4 (t(Fo) of 482 


Number of variables 


11 


RF = E\ J^a -Fa\/EFo 


0.0186 


Goodness of fit 


1.255 






Nd in 2a (0, 0, 0); 




Thermal displacement 


[A1 


Interatomic distances [A] 




Occupancy 


1.00(1) 


Nd: Utt = U22 = U33 


0.0482(5) 


Nd - 12 Sb 


3.4831 


Os in 8c (1/4, 1/4, 1/4); 




Thermal displacement 


[A1 


Interatomic distances [A] 




Occupancy 


1.00(1) 


Os: Uii = U22 = U33 


0.0025(1) 


Os - 6 Sb 


2.6239 


Sb in 24g (0, y, z); y: 


0.15597(3) 


Thermal displacement 


[A1 


Interatomic distances [A] 




z: 


0.34017(3) 


Sb: Uii 


0.0026(1) 


Sb - 1 Sb 


2.9033 


Occupancy 


1.00(1) 


U22 


0.0044(1) 


- 1 Sb 


2.9752 






U33 


0.0065(1) 


- 2 Os 


2.6239 










- 1 Nd 


3.4831 



fourth- and sixth-order terms of the angular momentum 
operators and W is an overall energy scale. The CEF 
contribution to the molar magnetic susceptibility can be 
determined from the expression^ 



Xcef(T) = iVA5>B 



E 



-2 E 




where Na is Avogadro's number, gj is the Lande g-factor, 



fiB is the Bohr magneton, pi 



-E,/{kBT) 



/Z is the ther- 



mal population probability (Z is the partition function), 
and the Ei's are the energies of the multiplets. However, 
the occurrence of ferromagnetic order in NdOs4Sbi2 re- 
quires the presence of a molecular field constant A to 
account for the effective field, with = Xcef ~ ^• 

The low-temperature value of fXcS (2.35 (Ib) indicates 
that the ground state of Nd^+ is either Fg or rg^-*. A 
wide range of xllw values with a Tq ground state can fit 
the x{T) data reasonably well (— 1 < a^LLW < 0.1). For a 
Fg ground state, a good fit to the x{T) data is only found 
for a;LLW ~ —0.4. All the fits imply that the splitting be- 
tween the ground and first excited states is greater than 
120 K. When compared with the p{T) data (to be dis- 
cussed later), the best fits are (I) xllw = —0.1, W — 3.9, 
and (II) xllw = —0.4, W = —7.2, corresponding to 
the following level schemes: (I) Fg (OK), F^^^ (180 K), 



F^^^ (420K)and(II)F^"^ (OK),F^'^ (220K), Fg (600K), 
with A = 1.39 Nd-mol/cm'^ and an effective ground state 
moment ~ 2.31/iB for both schemes. These fits are dis- 
played in Figs.^a) and (b). The Heisenberg interaction 
strength between a Nd ion and its eight nearest neighbors 
is estimated to be 0.522 K from A. 

Isothermal magnetization measurements M{H), dis- 
played in Fig.|2Ia), were made above and below the Curie 
temperature (6c)- In the paramagnetic state, the M (H) 
isotherms between 2 K and 5 K can be fit well by a Bril- 
louin function with effective J = 2.7 that includes a tem- 
perature shift of IK due to ferromagnetism. The M{H) 
isotherms in the vicinity of Qc were used to construct 
a conventional Arrott plot consisting of vs H/M 
isotherms, shown in Fig. |2tb) . The isotherms in the Ar- 
rott plot are linear and parallel in the vicinity of Qq for 
H < 1 T, indicating that NdOs4Sbi2 is a mean-field fer- 



.(2) 



(1) 




5 
4 

"e 

:^ 3 


E 

1 







4 5 




(c)/ 






- '2 inverse initial"/ 




from Arrott plot 









0.5 1 1.5 2 2.5 3 
H/M(T/u) 



Q. 



-1.76U ;f.u. 



(d) 



M -0.7(1 



S -4- 



2 
T(K) 



0.5 



1 1.5 
T(K) 



2.5 



Fig.2 Ho, P.-C.,etal, 

Weak ferromagnetism and heavy fermion behavior in single crystals of NdOs Sb 



FIG. 2: (a) Magnetization M vs magnetic field H isotherms 
and the saturation magnetization Mgat determined from M 
vs H isotherms below the Curie temperature Qc- (b) 
vs H/M isotherms (Arrott plot) for NdOs4Sbi2. (c) Inverse 
initial magnetic susceptibility (open circles) determined from 
= intercepts of vs H/M isotherms, compared to 
the Xdc^iT) data (pluses). The line is a Curie- Weiss fit. (d) 
H/M — intercepts of vs H/M isotherms plotted versus 
T. Positive values correspond to the spontaneous magnetiza- 
tion Man. 



romagnet. The intercepts of the linear fits to the {H/M) 
axis (= the inverse initial magnetic susceptibility) agree 
well with the low-temperature Xdc^(^) data (Fig. I^Jc)). 
The intercepts of the linear fit to the A'P axis are shown 
in Fig. I^d), where zero identifies 8c — 0.93 K. Below 
Qc, the intercept of the Af^-axis corresponds to the 
square of the spontaneous magnetization (Mg^p), which 
levels off and results in a small value of Mgp ~ O.T/ie/f-u.. 
However, a linear extrapolation of the negative M^-axis 
intercept back to zero temperature yields a much larger 
value of 1.76 /is/f-u., which is comparable to the satu- 
ration magnetization Mgat of 1.73 /ie/f-u. determined 
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Fig. 3 Ho, P.-C, et al.: Weak ferramagnetism and heavy fermion behavior ih single cryslals of NdOs ,,Sb^^ 

FIG. 3: (a) Low-temperature electrical resistivity p vs tem- 
perature T for magnetic fields H between OT and 6T for 
NdOs4Sbi2. The solid lines are guides to the eye. (b) Low- 
temperature electrical resistivity p vs T for magnetic fields 
H between 8T and 18 T for NdOs4Sbi2. The solid lines are 
guides to the eye. (c) Position of the shoulder Td in p{T) 
(which corresponds to Sc) vs H , with vertical bars indicat- 
ing the width of the transition at Td as described in the text, 
(d) High-T resistivity p vs T at H = OT, 0.5 T, IT, 2T, 4T, 
6T, 8T, 9T. 

directly from the M{H) data of NdOs4Sbi2 (Fig.EIa)). 
The value of Mgat is also consistent with the saturation 
magnetization found in NdFe4Pi2 (1.72/iB/f.u.)ii^ 



C. Electrical Resistivity 

Low-temperature electrical resistivity p{T) data for 
NdOs4Sbi2 in various magnetic fields from OT to 18 T 
are shown in Figs. |2Ia) and (b). The zero-field residual 
resistivity ratio RRR ee p(300 K) /p(0.02K) of ~ 45 indi- 
cates that the single crystal studied is of good metallurgi- 
cal quality (Fig. |3fa) and (d)). A shoulder occurs in the 
zero-field p{T) curve at ~ 1.2 K, below which p{T) has 
a sharp drop, indicating the development of an ordered 
state. The temperature Td at which this drop occurs is 
defined as the intercept of two lines, one of which is a 
linear fit to the data above the transition while the other 
is a linear fit to the data below the transition. The up- 
per and lower limits of the transition are defined as the 
temperatures midway between Td and the temperatures 
at which the data deviate from the linear fits. Shown in 
Fig. |3fc) is the field dependence of Td, which increases 
with increasing field up to 6 T, and is no longer observed 
in the p{T) data at higher fields. The temperature Td cor- 
relates with the onset of ferromagnetic order at 0.9 K de- 
termined from magnetization and specific heat measure- 




Fig.4 Ho, P.-C , etal. 

Weak ferramagnetism and heavy fermion behavior in single crystals of NdOs Sb 

FIG. 4: (a) Low-T electrical resistivity p vs temperature T 
with power-law (dashed line) and spin-wave (solid line) fits at 
various magnetic fields for NdOs4Sbi2. (b) Residual resistiv- 
ity po vs H. (c) Exponent n of the power-law fit vs H. (d) 
Energy gap Aspw from the spin- wave fit vs iJ". 

ments. Displayed in Fig. Efd) are the high-temperature 
p{T, OT < H < 9T) data for NdOs4Sbi2, which show a 
slight increase in p(T) with increasing field. 

In order to analyze the behavior of the resistivity be- 
low Qc, the p{T) data were fit with a power-law of the 
form p(T) = Po + BT". The fitting curves are plot- 
ted as dashed lines in Fig. 2{a) . The residual resistivity 
Po increases with increasing field and has a linear H- 
dependence above 8 T (Fig.lJb)). Between OT and 18 T, 
the exponent n varies from 3 to 4, which indicates that 
NdOs4Sbi2 exhibits neither typical Fermi-liquid (n ~ 2) 
nor typical non- Fermi-liquid (n < 2) behavior (Fig.^fc)). 
Since ferromagnetic ordering occurs below Oc, electron- 
spin wave scattering was considered with the form^ 

p{T) = PO + A-^ (l + 2-f-) exp (~^) , (3) 

where Agp^ is the spin wave energy gap, which may re- 
sult either from magnetic anisotropy or from broken sym- 
metry due to presence of a CEF. This formula describes 
the p(r, H) data well, and the fitting curves are shown as 
solid lines in Fig.^a). As determined from these fits, the 
spin- wave energy gap Agp^ is ^ 0.75 K at OT, increases 
approximately linearly to ^ 4.5 K as the field increases 
to 12 T, and then drops to - 3.4 K at 18 T (Fig. g^d)). 

Figure l^a) displays the zero-field p(T) data, which 
have a slight negative curvature at ~ 130 K that may 
be related to scattering from the CEF levels. In order 
to analyze the CEF contribution to p(T) at OT, it is 
necessary to subtract a lattice contribution piat(7') and 
an impurity contribution Pimp 9.4 pfl-cm) from the 
resistivity data. Usually, pia.t{T) is estimated from an 
isostructural nonmagnetic reference compound; in the 
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FIG. 5: (a) Zero-field resistivity p and the estimated piat+pimp 
vs temperature T for NdOs4Sbi2, where pimp ~ 9.4 [£l-cm. 
(b) Temperature dependence of the incremental resistivity 
Ap = p — Plat — Pimp and CEF fits for two different ground 
states: (dashed line) and Fg'^'' (solid line), where a;/(l — x) 
is the ratio of s-f exchange to aspherical Coulomb scattering 
for NdOs4Sbi2. 



case of NdOs4Sbi2, the resistivity of LaOs4Sbi2, which 
has an empty 4-f shell, was used. However, above 100 K, 
/9(T) of LaOs4Sbi2 exhibits a significant negative cur- 
vature that is common in La-based compounds such as 
LaAl2i^ii2S This curvature is generally less pronounced 
in Sc-, Y-, and Lu-based compounds, which have com- 
pletely empty or filled 4f-electron shells. However, the 
compounds ScOs4Sbi2, YOs4Sbi2 and LuOs4Sbi2 have 
not yet been synthesized. Above 100 K, piat(r) was es- 
timated from SmOs4Sbi2, as p(T) of SmOs4Sbi2 has an 
approximately linear-T dependence between 100 K and 
300 KiS The s-f exchange scattering effect in /9(T) of 
SmOs4Sbi2 only appears below 80 K due to the small en- 
ergy splitting (~ 35 K) between the ground and the first 
excited states in that compound. Thus it is reasonable to 
assume that p{T) of SmOs4Sbi2 for lOOK < T < 300 K 
is due to electron-phonon scattering and use it as an 
estimate of the high-temperature portion of piat(T') for 
NdOs4Sbi2. 

The incremental resistivity Ap(r) = p(r) — piat(r) — 
Pimp (Fig. Elb)) is best described by two energy-level 
schemes that are consistent with the CEF analysis of 
x(r) discussed previously. During the CEF analysis 
of Ap(T), it was also found that s-f exchange scatter- 
ing alone could not entirely account for Ayo(T); other- 
wise. Pimp would always be negative, which is unphysical. 
Thus, the effect of aspherical Coulomb scattering2ii2Si2& 
was also considered, with the ratio between the s-f ex- 
change scattering and the aspherical Coulomb scattering 
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Weak ferromagnetism and heavy fermion behavior in single crystais of NdOs^Sb^^ 

FIG. 6: (a) Isotherms of electrical resistivity p vs magnetic 
field H for NdOs4Sbi2. (b) and (c) Calculated p ys H 
isotherms including the effects of s-f exchange and aspheri- 
cal Coulomb scattering using the CEF parameters determined 
from zero-field data. 



defined as a; : (l-x). In Fig. Elb), the fit curves represent- 
ing the two best-fit energy schemes are plotted in com- 
parison with Ap{T), which is described well by both fits 
above 40 K. The departure of Ap(T) from the fits below 
40 K may be due to a reduction of the electron scattering, 
resulting from the development of the coherent heavy- 
fermion ground state in NdOs4Sbi2 as the temperature 
is lowered. The p{H) isotherms are displayed in Fig. Efa) 
and qualitatively agree with the p{H) isotherms gener- 
ated using the cubic-CEF parameters determined from 
the zero-field fits to Ap(T) (Figs.l^b) and (c)). 

Measurements of p{T) were performed from 1 K to 300 
K under nearly hydrostatic pressure P between kbar 
and 28 kbar (Fig. [7|l . In Fig. [7| it can clearly be seen 
that the pressure-induced change in the high-T electrical 
resistivity is much more pronounced than the resistivity 
change induced by high field (Fig.|3{c)), in the pressure 
and temperature ranges of this investigation. However, 
at low temperatures, the value of is more strongly 
influenced by an increase in fleld H (Fig. Etb)) than by 
variation of P (Fig. [Tlinset (b)). 



D. Specific Heat 

Specific heat divided by temperature C/T vs T data 
for NdOs4Sbi2 are shown in Fig.lSJa). The data reveal a 
pronounced peak at ^ 0.8 K that correlates well with 
the magnetic ordering temperature Qq inferred from 
the shoulder of p(T), the divergence of Xdc{T), and the 
Arrott-plot analysis. No obvious Schottky anomaly as- 
sociated with CEF energy splitting was observed in the 
C(T) data below 20 K. This is consistent with the fact 
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FIG. 7: Electrical resistivity p vs temperature T for 
NdOs4Sbi2 at various pressures P up to 28kbar. Inset (a): 
Expanded view of the resistive ferromagnetic transitions. In- 
set (b): Temperature of the drop in p{T) due to the onset 
of ferromagnetism Td (approximately corresponding to the 
Curie temperature Be) vs P with vertical bars indicating the 
width of the transition. 



that all the fits to the dc magnetic susceptibility data 
imply that the CEF splitting between the ground and 
first excited states is greater than 120 K. The Schottky 
anomaly due to 120 K CEF splitting would exhibit a 
peak at ~ 40K and make a negligible contribution to 
C{T) 20 K. Such a peak would also be difficult to resolve 
against the large lattice background. 

In Fig. IHa), the specific heat of LaOs4Sbi2 is dis- 
played in comparison with that of NdOs4Sbi2, where 
the electronic specific heat coefficient 7 and the Debye 
temperature Od of LaOs4Sbi2 are ^ 36mJ/mol-K^ and 
^ 280 K, respectively. After the specific heat of nonmag- 
netic LaOs4Sbi2 (an estimate of the lattice heat capac- 
ity of NdOs4Sbi2) is subtracted from the specific heat 
of NdOs4Sbi2, the difference divided by temperature 
5C /T is plotted against and shown in the inset of 
Fig. IS^a). The value of 7 for NdOs4Sbi2 estimated from 
the 6C/T vs plot ranges from ~ 436mJ/mol-K^ to 
- 530mJ/mol-K2. Below 13 K, 5C/T vs is not con- 
stant, which could be due to a difference between the ac- 
tual lattice heat capacities of NdOs4Sbi2 and LaOs4Sbi2. 
Nevertheless, the curvature in C/T vs of NdOs4Sbi2 
(inset of Fig.|Sfb)) and the magnetic transition occurring 
at ~ 1 K cause difficulties in the analysis of the data us- 
ing the typical formula C/T = 7-1- (iT"^ . Since we have 
strong evidence against the possibility of a CEF Schot- 
tky contribution from the analysis of the x(T) data, the 
curvature in C/T vs T^ is possibly due to either the rat- 
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Weak ferromagnetism and heavy fermion behavior in single crystals of NdOs Sb 

FIG. 8: (a) Zero-field C vs T for NdOs4Sbi2 and LaOs4Sbi2 
below 20 K. Inset: 5C /T vs below 20 K, where SC = 
C(NdOs4Sbi2) - C(LaOs4Sbi2). The extrapolated values of 
the two dashed lines at OK set the lower and upper limits 
of 7(NdOs4Sbi2) - 7(LaOs4Sbi2). (b) A comparison between 
C of NdOs4Sbi2 and a fit of + Ciat, where Ca = 7T 
and Ciat is composed of Cneb and Ceiii as described in the 
text. The electronic specific heat coefficient 7, the Debye 
temperature (6d), the Einstein temperature (0E),and the 
coupling constant r for the Einstein oscillator are estimated 
as 520mJ/mol-K^ 255 K, 39K and 0.48 respectively. Inset: 
C/T vs for NdOs4Sbi2 below 20 K for NdOs4Sbi2. 



tling motion of the Nd atoms or a narrow peak in the 
density of electronic states near Fermi energy, such as 
a Kondo resonance^ However, we do not consider the 
application of the resonance level model (RLM)^ to be 
appropriate, because the magnetization data above Qq 
are fit well by a Brillouin function and no obvious Kondo 
effect is observed in the p{T) data of NdOs4Sbi2. In 
the previous specific heat studies of the Tlo.22Co4Sbi2 
filled skutterudite compound, Sales et aljS found that 
the difference in heat capacity between Tlo.22Co4Sbi2 
and the unfilled skutterudite compound Co4Sbi2 can be 
accurately described by a quantized oscillator (Einstein 
model) with an Einstein temperature Oe of 55 K. Since 
the X-ray structural refinement at room temperature for 
NdOs4Sbi2 indicates a small 9e 45K associated with 
the Nd atoms, it can be assumed that the Nd atoms par- 
tially act like Einstein oscillators with a mixing ratio r, 
and the lattice contribution to the specific heat can be 
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Weak ferromagnetism and heavy fermion behavior in singie crystais of NdOs Sb 



FIG. 9: (a) Incremental specific heat AC (AC = C — Cd — 
Ciat) (left axis) and the magnetic entropy Smag (right axis) 
vs temperature T. (b) Logarithmic plot of the power-law fit 
(dotted line) and the anisotropic-spin-wave fit (dashed line) to 
the incremental specific heat AC(T) after the electronic and 
lattice contributions have been removed (in a very limited 
measuring range below 0c)- 



expressed as Ciat = CEin + Cuch, 

(eE/r)2e(«'^/^) 



Ceii 



3R- 



[c(eE/T) _ 1]5 



Cocb = (17 - r) 



12^4 / T 



(4) 



(5) 



where R is the universal gas constant, 8d represents 
the Debye temperature, and r < 1. The effective De- 
bye contribution is (17 r)/r times bigger than that of 
the Einstein-like Nd rattling motion due to the partici- 
pation of the rest of the atoms in the unit cell. Below 
20 K, a least squares fit of Cd + Cut to the C(T) data 
was performed, where Cei — 7T is the electronic specific 
heat, from which estimated values of 7, 0d,0Ej and r 
were derived. The fitting curve is plotted in Fig. |Hl[b) 
along with the C(T) data of NdOs4Sbi2 as a compari- 
son. The values obtained for 7, 6d, Be, and r are es- 
timated as 520mJ/mol-K2, 255 K, 39 K, and 0.48, re- 
spectively. The value of 8d is close to the Debye tem- 
perature of LaOs4Sbi2, the value of 8e is comparable 
to that determined from the X-ray data, and the value 
of 7 is close to that estimated from the simple subtrac- 
tion of the LaOs4Sbi2 specific heat data, suggesting that 
NdOs4Sbi2 is possibly a heavy fermion compound. 



The temperature dependence of the magnetic entropy 
'S'mag was derived from the integration of AC/T vs T and 
is shown in Fig. EI a) , where AC = C — Cd — Ciat ■ The 
magnetic entropy (S'mag) in NdOs4Sbi2 reaches 0.69R (« 
Rln2) at G.85K and levels off at a value of ^ 1.14R 
(« Rln3). The magnetic entropy reaches ^ 74% of its 
full value at Oc, and a noticeable magnetic contribu- 
tion persists up to ~ 3K, revealing the existence of mag- 
netic fluctuations above 8c- It has been argued for 
the antiferromagnetic system Gdi_3;Y2.Ni2Si2 that mag- 
netic fluctuations can contribute to C(T) at tempera- 
tures up to 5 times the Neel temperature)^ Thus, we 
cannot completely rule out the possibility that the short 
range magnetic correlations near the ferromagnetic tran- 
sition temperature regime give rise to enhancement of 
the specific heat of NdOs4Sbi2i^S*2Ma However, it IS un- 
likely that they would account for a large fraction of 
the enhanced specific heat, due to the following argu- 
ments: (i) The paramagnetic-state M[H) isotherm data 
(2 - 5K) scale well with a Brillouin function. The M{H) 
data at 2K and 3K, and the Xdli"^) data along with 
the initial from the Arrott Plot analysis do not 
show obvious evidence of magnetic fluctuations above 
2K (displayed in Fig. [3 (a) and (b)). (ii) The value 
estimated for 7 (520 mJ/mol-K^) is within the upper 
and lower limits of 7 (436 - 53GmJ/mol-K^) estimated 
from the analysis done by comparing the specific heat 
of NdOs4Sbi2 with that of the nonmagnetic compound 
LaOs4Sbi2 suggesting that our analysis is justified, (iii) 
The lower limit (4K) of the fitting range used to deter- 
mine 7 from the formula C — Cd + Ciat is four times the 
Curie temperature, which is safely higher than the tem- 
perature at which the calculated Smag saturates, (iv) 
Heavy fermion behavior has been found in the neighbor- 
ing compounds PrOs4Sbi2 (7 ~ 600 mJ/mol-K^)'*^ and 
SmOs4Sbi2 (7 « 880 mJ/mol-K^)f2i earlier studies of 
the related compound NdFe4Sbi2 also reported possible 
evidence for an enhanced electron mass.^^'^'^ Consider- 
ing all of these points, it seems likely that NdOs4Sbi2 
displays heavy fermion behavior. 

Since S'mag lies between Rln2 {— 0.69R) and Rln4 {— 
1.39R), it is difficult to determine conclusively whether 
the Fg doublet or Fg^^ quartet is the Nd'^^ ground state. 
If the Fg doublet is the ground state, then the extra en- 
tropy may result from another degree of freedom, such as 
a tunnelling mode or off-center mode of Nd^+ ion in an 
Sb-icosahedron cage.^^^ If the Fg^-* quartet is the ground 
state, the smaller entropy may be due to an overesti- 
mated lattice contribution to the specific heat or transfer 
of entropy to the conduction electron system. 

Figure Efb) displays the incremental specific heat 
AC(T) after the electronic and lattice contributions have 
been removed. Even though the range of the AC(r) data 
below 8c is very limited, the data were fit with spin- 
wave formulas AC(T) oc T" for magnetically isotropic 
metals and AC(r) oc T^/"^ exp(-Aspw/T) for magneti- 
cally anisotropic metals.^* From the ffist formula, the 
value of the exponent n « 2.3 is higher than the value of 
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1.5 expected from a spin wave in an isotropic metal. The 
spin-wave energy gap Agp^ determined from the second 
formula is ^ 0.54 K, consistent with the value of 0.75 K 
determined from the zero-field resistivity data. 

IV. SUMMARY 

We have performed X-ray diffraction, electrical re- 
sistivity, magnetization, and specific heat studies of 
NdOs4Sbi2 single crystals, which exhibit interesting 
strongly-correlated-electron behavior. X-ray experi- 
ments have revealed full occupancy of Nd sites and a 
large mean square displacement of Nd ions in NdOs4Sbi2. 
The compound NdOs4Sbi2 exhibits mean-ficld-type fer- 
romagnetism with Oc ~ 0.9 K. The value of 7 es- 
timated from the analysis of the specific heat is 
large, ~ 520mJ/mol-K^ {m* ~98me), indicating that 
NdOs4Sbi2 is a possible heavy fermion compound. A 
cubic CEF analysis suggests two best-fit energy level 

schemes: (I) re (OK), v'^^ (180 K), (420 K): (II) 
(OK), (220 K), Te (600 K), both with a molecu- 
lar field parameter A = 1.39 Nd-mol/cm^. The electrical 
resistivity data indicate that both s-f exchange and as- 
pherical Coulomb scattering are present in NdOs4Sbi2. 



Low-temperature electrical resistivity data suggest the 
possible existence of spin- wave excitations below Oq- 
The uncertainty in the CEF-energy-level scheme ground 
state and the possible existence of spin-wave excitations 
may be resolved by further neutron scattering experi- 
ments. 
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